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Q-t I. INTRODUCTION 

Since the early days of Quantum Mechanics the problem of finding potentials such that the one-dimensional 
Schrodinger equation (ODSE) 

^ + [ fe 2 _ y(x) ] u{x) = 0) k 2 = ^, xe (-00,00) (1) 

admits exact solutions attracted the attention of several physicists and mathematicians 1]. The importance of the 
problem is directly related to the fact that an exact knowledge of the energy spectrum and of the wave functions 
of the ODSE turns out to be very useful in many applied problems @, Q ■ It is well-known that the ODSE can be 
solved analytically for a large number of potentials provided that we are able to hnd a coordinate transformation 
y = y{x) and eventually a change of the dependent variable u such that ([1} can be transformed into a linear ordinary 
differential equation of mathematical physics. Here, we begin by reviewing the problem of reducing the ODSE to the 
^vq \ hypergeometric equation 



y^-y)— + [c~(a + b+l)y} — -abv(y) = 0, yelcR (2) 



dv 

' dy 2 ' L v ' ,ai dy 
witha,6,ceC. \M 

considered only the case when the potentials depend explicitly on the variable x. This scenario 
has been further extended by [8, 9], where more general classes of potentials exhibiting an implicit dependence on 
the coordinate x have been derived. Note that the so-called Ginocchio potentials [l(| are a special case of the more 
general class of solvable potentials obtained by [9(]. Moreover, derived a large class of potentials which include 
the Natanzon class by considering the case when a second order linear differential equation can be connected to a 
ODSE by means of Liouville transformations. We emphasize that the work of does not cover the case of the Heun 
equation, since the author assumed that the Bose invariant is a rational function with numerator represented by a 
polynomial of degree two or less. As we will see in Section ITVl such a numerator becomes a polynomial of degree four 
in the case of the Heun equation. Our main contribution is the derivation of the most general potential such that the 



ODSE can be transformed into the Heun equation [12 1 

d v f 7 (5 e \ dv a/3y 



d V 2 \V y-l y-aj dy y(y-l)(y-a) 



-v(y) = 0, y € I C R, (3) 
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where a £ C\{0, 1}, e = a + /3 + l — 7 — S, and a, (3, 7, S, and q are arbitrary complex parameters. For the role 
played by in General Relativity we refer to [li|. Since the hypergeometric equation is a special case of the Hcun 
equation with q = a(3a and e = 0, our potential generalizes the Natanzon class. We also treat the same problem for 
the confluent, biconfluent and triconfluent Heun equation. Moreover, we also find the most general class of potentials 
such that the ODSE can be transformed into the generalized Heun equation (GHE) [14| 



dy- 



+ 



1 



fJ-o 



1 



/'1 



1 



I) 



y-a 



dv f3 + p lV + fcy 2 
dy y(y-l)(y-a) 



v(y) = 0, 



(4) 



where a G C\{0, 1} and no, fj,i, fi 2 , a, /3 Q , j3i, /3 2 are arbitrary complex numbers. To underline the importance of 
equation (0} we recall that it contains the ellipsoidal wave equation as well as the Heun equation (a — f3 2 — 0) and 
thus the Mathieu, spheroidal, Lame, Whittaker-Hill and Ince equations as special cases. For physical applications 
of (j4} we refer to [la ]. Our work is organized as follows. In Section |H] we outline the derivation of the potential 
found by @ , extend the findings thereof by computing two new classes of potentials and derive a non trivial condition 
under which the potential is symmetric. In Section IlIII we shortly review the derivation of Natanzon's potential and 
show that all Iwata's potentials are contained in the Natanazon class. In Section HVl we compute the most general 
potential such that ODSE can be reduced to the Heun equation. This potential will be called potential of the Heun 
class. Starting from this result we also derive expressions for the potentials when the ODSE can be transformed into 
a confluent, doubly confluent, biconfluent and triconfluent Heun equation. We conclude our work by deriving the 
most general potential such that the ODSE can be transformed into the generalized Heun equation. Since the Hcun 
equation is a special case of the GHE, we obtain a further generalization of the potentials belonging to the Heun class. 
A comment on the choice of parameters entering in the potentials derived in the next sections is in order. From here 
on, we will be concerned with real parameters in agreement with the choice of a real potential in the ODSE. However, 
the choice of complex parameters would also be possible in the ODSE when the potential is complex, a case which 
has some relevance in the physical literature and goes under the name of optical potential. 

II. IWATA'S POTENTIAL 

We start as in Q by restricting our attention to those potentials V such that (TTJ) can be reduced to ^ after a 
suitable transformation of the coordinate x and of the dependent variable u. If we transform ([2]) according to 



v = y{x), v(y) = s(x)u(x), 



we obtain 
dx 2 



s y c- (0 + 6+1)3/ / 
2 r + 7Z r V 

s y y(i-y 



du, 
dx 



s y s c — (a + 6 + l)y s > ab(y ) 2 
s y s 2/(1 - y) s V y(l- y) 



In the process of reducing © to (fTJ) we must require that 



S y . c - a + b + l)y , 

2 rH t, \ y =0. 

s y y{±-y) 



(5) 



i(x) = 0. (6) 



(7) 



At this point (JT)) can be used to eliminate the function s in the last bracket in (J6j) and we end up with the equation 



where 



with 



— -7: + U(x)u(x) = 0, 
dx z 



U(x) = T(x) + A(a, b)X(x) + B(a, b, c)Y(x) + C(c)Z(x) 



(8) 



Y(x) = 



' \2 



Z(x) = 



(y) 



'\2 



and 



4j/(l-2/) 2 ' w 4y 2 (l-2/) 2 ' 
A(a,b) = 1- (a-b) 2 , B(a,b,c)=2[c(a+b-l)-2ab], C(c) = 2c-c 2 . 



(9) 



(10) 
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Since the parameters a, b and c may depend on the quantity k defined in fT}, we can assume without loss of generality 
that 

U(x) = T{x) + A(k)X(x) + B(k)Y(x) + C{k)Z{x). (11) 

Differentiating the above expression first with respect to k and then with respect to x, we end up with the auxiliary 
equation 

Ax' + BY' + CZ = 0, (12) 
where ' = d/dk and = d/dx. We must distinguish among the following three cases 

1. The functions X , Y and Z are linearly independent. Then, A = B = C = meaning that A, B and C are 
numerical constants. This case is not of practical relevance. 

2. Only one function is linearly independent. For instance, there exist constants, say p and a such that Y = pX 
and Z = crX . A simple integration gives Y = pX + c\ and Z — aX + C2 . On the other side by means of © 
we find that the coordinate transformation y should satisfy at the same time the differential equations 

Af , 1 ( ~ ~ p) (y ) 2 = c i> 77T~^\ i ~2 ~ a ) (y ) 2 = C2 ' 
4(i -r) \y J 4(i -y 2 ) vr / 

This is the case if and only if y — 0, but then y is a constant function. Hence, this case must be disregarded. 

3. Two of the functions X , Y and Z are linearly independent. We have three distinct cases. 
Let Y and Z be linearly independent. Then, there exist constants p and a such that 

X' + pY' + aZ' = 0. (13) 

By means of ([13]) we can rewrite fTS]) as follows 

(B - pA)Y' +{C- oA)Z = 0. 

Since Y and Z are linearly independent, it must be B — pA — and C — a A = 0. Integrating these equations with 
respect to k we obtain 

B = pA + pi, C = aA + ax, pi, ai el (14) 

and U can be expressed as follows 

U(x) = T(x) + A [X(x) + pY(x) + aZ{x)] + p x Y{x) + a x Z{x). (15) 



If we integrate (j 1 3(1 with respect to x and make use of (|9|), we discover that the coordinate transformation y must 
satisfy the following first order nonlinear autonomous differential equation, namely 

j ■" ! ' 1 ~K, K G M\{0}, R(y) = y 2 + py + a. (16) 



4y 2 (l-y) 2 

Taking into account that U = k 2 — Vj and applying (| 15[) . we find that 



V I (x)=k 2 +(jL)j -A[X(x)+pY(x)+aZ(x)]-p 1 Y(x)-a 1 Z{x). 

With the help of © and (jTS]) it is straightforward to check that 

A [X(x) + pY(x) + oZ(x)\ + Pl Y(x) + u x Z{x) = k fl - (o - b) 2 } + K (^ + gl ) . 

R{y) 

Since the potential cannot depend on k, the parameters a and b must satisfy the following condition 

[a-b) 2 = - k - (17) 
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and the expression for the potential reduces to 



Vj(x) 



y 

2y' 



v" \ _ y 2 + (p + pi)y + <j + <7i 
2y' K R(y) 



The first two terms on the r.h.s. of the previous expression can be computed as follows. First of all, let 

W(l-2/) 2 



g(y) 



R(y) 



(18) 



Then, (1161) can be written in the compact form (y ) 2 = g(y) from which we can easily derive the useful identities 



Taking into account that 



" = ldg y ^ld 2 g 
V 2dy y 2dy 2 ' 



3 y 



we conclude that 



y 

V 



Vl{x) = W 9 {Ty) 



ly 3(/) 2 ly" 



V / 4 I y J 2 y Ag 2 y 



3 ( dg\ ld 2 g y 2 + (p + p x )y + a + a x 



R(y) 



with y = y(x) and g given by (|18[) . If we replace (|18|) in the above expression, we obtain as in 



Vj{x) 



ny 2 (l - yf 

R{y) 



(l-2y)(2y + p) 5(2y + pf 



y 2 (l-y) 2 R(y) y(\-y)R(y) 4R 2 (y) 



y 2 + (p + pi)y + o- + o\ 

R(y) 



(19) 



We call Vi Iwata's potential in order to distinguish it from a similar result obtained in Note that (fT9")) reduces to 
the modified Poschl- Teller potential barrier [29 



V(x) 



V 



cosh (ax) 



a>0 



(20) 



whenever 



i n 2 ^3 

p=-l, o = 0, K = a pi = cti = -pi + -. 

a z 4 

The potential V% will vanish identically whenever 

p = -2, cr = 1, (01=0-1= 0, 



p = -l, ct = 0, pi = 0, o-i = -, 

or 

p = (7 = pi = (7l = 0. 

Relations connecting the parameters of the hypergeometric equation with those appearing in ()19[) can be derived as 
follows. Since the energy of the particle must be real, from (fTT|) we have (a — b) 2 < if k > 0, and (a — o) 2 > if 
k < 0. Thus, for k > we obtain 



a — o = 




5 



which implies that the real parts of a and b coincide. Moreover, for re < we get 

k 



a — b 



implying that the imaginary parts of a and b coincide. From (flU)) and ([T4"f we have 



2[c{a + b- 1) - 2ab] = p[l - (a - b) 2 ] + p 1: 
2c- c 2 = a[l - (a - b) 2 } + ai. 



Replacing ([T7| into ([23]) we get 



fc 2 

c = 1 + \/ 1 — (7 — (Ti — a — . 



(22) 
(23) 



(24) 



Taking into account that 



can be manipulated as in [8[ to give 
(1 + cr + p)(a - b) 2 - (a + b - c) 2 = p + pi + cr + ai, 



we find that 



6 = 2 



1 + \j 1 — cr — ai — cr— + % | y/J + Pi+c + ci 



1 +P + CT 



fc 2 



1 — cr — (7i - (7— + i ( \/ /) + /), — rr -!- rr. 



1 +/5 + CT 



fc 2 - — 
\/K 



if re > 0. In the case re < we get 



1 + 



k 2 . l + p + a J9 

1 - cr — cri - (7 h i\ p + pi + a + ax -\ k z 

rev K 



»=2 1 



k 2 . l + p + a l0 

1 — (7 — (7l — (7 h «W p + pi + <J + G\ H fc Z 

rev re 



(25) 
(26) 

(27) 
(28) 



We point out that Q treated only the case when Y and Z are linearly independent. In what follows we extend the 
work of 1 by considering two additional cases generating new classes of potentials. To this purpose, let X and Y 
be linearly independent. Then, there exist constants p and a such that Z + pX + uY = 0. This condition together 
with (fT2")) implies that A — pC = and B — aC = 0. Integrating these equations with respect to k we get A — pC + p2 
and B = crC + 02 with p2,&2 € K. At this point U can be expressed as follows 



U{x) = T(x) + C [pX(x) + dY(x) + Z{x)) + p 2 X(x) + a 2 Y(x). 



(29) 



If we integrate (|13[) with respect to x and make use of @, we discover that the coordinate transformation y must 
satisfy the following first order nonlinear autonomous differential equation, namely 



(y') 2 R(y) 



= K, reeM\{0}, R(y) = py 2 + ay + l. 



4y 2 {l-y) 2 

Taking into account that U — k 2 — Vu and applying (|15p . we find that 

/ \ 2 / „ \ ' 



(30) 



V n {x) 



V 



V 



-hrr -C[pX{x) + aY{x) + Z{x)]-p 2 X(x)-a 2 Y{x). 



2y 



22/ 



With the help of (JQJ) and ([30]) it is straightforward to check that 



C [pX(x) + oY(x) + Z{x)} + P2 X{x) + a 2 Y(x) = re(2c - c 2 ) + 



K{p 2 y 2 + <J 2 y) 
R(y) 
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Since the potential cannot depend on k, the parameter c must satisfy the following condition 

k 2 

2c-c 2 = ^ (31) 

K 

and the expression for the potential reduces to 

' it \ 2 / " \ ' ~ 

y \ v \ K (p2y 2 + a iy) 



V n (x) 



2y' J \2y' R( y ) 



The first two terms on the r.h.s. of the above expression can be computed with the same procedure we adopted for 
the case when Y and Z are linearly independent. Hence, we end up with a new potential 



V u (x) = 

R(y) 



1 2p (l-2y)(2py + a) 5{2py + a) 2 



y 2 (l-y) 2 R(y) y(l-y)R(y) 4R 2 (y) 



~ K P < +CJ2y . (32) 

R(y) 



The parameters a, b, and c can be expressed in terms of the parameters of the potential and k as follows 



c + O 1 / _fc 2 k 2 / Jk 2 

a=—— h -\ 1 - p 2 - o 2 - (1 + p + cr) — , b = a-fl, c=l + \l-—, fi = \1 - p 2 - p— . 

22y k Vkv k 

It is interesting to observe that the potentials Vj and Vjj will coincide if and only if 

/J = l=/92, V = P, K = K, <72 = P + Pl, cr = 1, O"! = — 1. (33) 

Moreover, taking into account that in the case of the modified Poschl- Teller potential the coordinate transformation 
y is given by y(x) — cosh 2 (ax) with a > 2], it is not difficult to verify with the help of (|30|) that this potential does 
not belong to the class Vj/. A further class of potentials, here denoted as Vm can be obtained by considering the case 
when X and Z are linearly independent. Then, there exist constants p and a such that Y + pX + aZ = 0. If we 
use this condition together with (TT2"j) . we find that A — pB = and C — aB = 0. After integration of these equations 
with respect to k we get A = pB + p^ and C = aB + 03 with P3, 173 G K. Also in this case U can be expressed as 

U(x) = T{x) + B [pX(x) + Y(x) + aZ(x)\ + p 3 X(x) + <r 3 Z(x). (34) 

If we integrate (1131) with respect to x and make use of ©, we discover that the coordinate transformation y must 
satisfy the following first order nonlinear autonomous differential equation, namely 

4y™y)2 R(y) = py 2 + y + °- (35) 

Since U — k 2 — Vj/j and applying dT5|) . we find that 

Kj 7 (a:) = A; 2 + ~ (|/) ~ B & X W +dZ W ~ P* X W ~ a3Z ( x ^ 

With the help of © and (|55|) it can be checked that 

B \pX(x) + Y(x) + aZ(x)} + p 3 X(x) + a 3 Z{x) = 2k [c{a + b - 1) - 2ab] + ^1 + (J3 \ 

R{y) 

If we choose the parameters a, b and c so that 

k 2 

c(a + b-l)-2ab= — , (36) 
2k 

we ensure that the potential Vm will not depend on the energy of the particle. Hence, we obtain 

VlIl[x) = ll\\ll\ W + r^ 

W W R{y) 



7 



Computing the first two terms on the r.h.s. of the above expression as we did in the previous cases, we end up with 
the following potential 



T/ / ^ % 2 (l - y) 2 



1 2p , (l-2y)(2py + l) 5(2py+l) 5 



y 2 (l-?/) 2 R(y) y(l-y)R(y) AR 2 (y) 



^ + fT3 . (37) 

R(y) 



The parameters a, 6, and c can be expressed in terms of the parameters of the potential and k 2 as follows 



c + fi 1 / fc 2 ~ / k 2 ~ k 2 

a = — — i- 2Y 1 - ^3 - ^3 - (i + p + ^Y' & = a_fi > c = 1 + Y 1 - CT 3 - n= yi- P3~p-~- 

Already comparing the last two terms in (l3"2"j) and (|37p we see that in general the potentials Vn and Vm do not in 
general coincide. They will coincide if and only if 



K = K, p = p, (7 = 1= (7, (72 = = (73, p2 = P3- (38) 

Moreover, the potentials Vj and Vj/j will coincide whenever 

K = K, p = p = p 3 = l, p 1 = -l, a = a, a + a 1 =cr 3 . (39) 



Finally, by means of (|33|) , (|38[) and (1301) it can be easily verified that there is only one potential living in all three 
classes. This happens when 

K — K = K, p = p = p = p2 = p3 = l, pi = <7i = — 1, (7 = (7 = (7 = 1, (7 2 = (7 3 = 0. 

Trivially, Iwata's potential becomes an even function whenever the coordinate transformation y is even. Is it possible 
for Vi to be symmetric even though the coordinate transformation is not? The next result shows that if y is non- 
symmetric and a certain set of functions is linearly independent, the potential must be zero, i.e. the only symmetric 
potential for such a choice of y is the trivial potential. 

Theorem II. 1 Let D C R be a common domain of definition for y and v with v(x) — y{—x) and suppose that 
y(x) 7^ v{x) for at least one x G D C M. If the functions in the set IA = Ui=i with U\ = {y n — v n \ n = 1, • • • , 6}, 
U-2 = {v 6 y n — y 6 v n | n = 1, • • • , 5} 7 and = {y m v n — v m y n \ m, n = 1, • • • , 5; m > n} are linearly independent, then 
the potential Vi as given in \19\) vanishes under the assumption that it is symmetric, i.e. Vi(x) = Vi{— x). 

Proof. If we rewrite the potential as 

Vl{x) = -TTY P(v) = jl^y n , q(y) = R 3 (y) = Y / b ny n 

n=0 n=0 

with 

ao = 4ct 2 (ct + (7i — 1), a\ = 4a(pia + 3pa + 1o\p — 3p), 

a 2 = 12p(7 + 12(7 2 + 8p/9l(7 - 3,0 2 + 4(7l/9 2 + 12p 2 (7 - 24(7 + 8(71(7, 

a 3 = 16/?cr + 40(7 + 8(7ip + 4pip 2 + 8pi(7 + 4p 3 + 2p 2 -8p, 
a 4 = 9p 2 + 4cti - 12cr + 8pip + 12p, a 5 = 4p lr 

bo = ct 3 , bi—3pcr 2 , 62 = 3p 2 (7 + 3cr 2 , b 3 — p 3 + 6per, &4 = 3cr + 3p 2 , 
65 = 3p, b 6 = l 

and set v = y(—x), then Vj will be symmetric, i.e. Vi(x) = Vi(— x) if and only if 

p{y)q(v) -p(v)q(y) = 0. 

This relation can be equivalently written as 

5 55 
Y,(b a n - a b n )(y n - v n ) - a (y 6 - v 6 ) + ^ a n (v 6 y n - y 6 v n ) + ^ a m b n <p mn = 0, p mn = y m v n - v m y n , (40) 

n— n— 1 m.n— 1 
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where we made use of the fact that be = 1. Since the last term in the above equation vanishes for n = m and 
fmn = -^ntii! we finally obtain the following more economic version of (j40[) 

5 5 5 

J2( b o a n - a b n ){y n -v n )-a (y 6 -v 6 ) + Y,an(v 6 y n -y 6 v n )+ ^ a m b n <p mn = 0. (41) 

n— n— 1 m,n— 1 



Let us introduce the sets of functions U\ — {y n — v n \ n = 1, ••• ,6}, Ui — {v 6 y n — y 6 v n \ n = 1, ■ • • ,5}, and 
W3 = {y TO u n — v m y n I m, n = 1, • • • , 5; m > n}. If the functions in the set U = [j i=1 U% are linearly independent, we 
can conclude that (141j) will hold true for every x whenever all coefficients in (|4"Tj) vanish, i.e. 

a = 0, 6 a„ - a &„ = V n = 1, • • • , 5, a\ = ■ ■ ■ = a 5 = 0. (42) 

Note that the last condition implies that a m b n = for any m,n = 1, • • • ,5. Since do = ai = • • • = (25 = 0, then 
= and thus Vj must vanish. □ 



Remark II. 2 It is obvious that U6\) has a stationary solution y\ = 0. Hence, among all possible nonsymmetric 
coordinate transformations those such that y(—x) = —y{x) must be ruled out since they would intersect the solution 
yi and this would violate the Existence and Uniqueness Theorem for first order differential equations. 

The system (|4"2"| has been solved with Maple 15 and as a result we found that the only solution is given by p = a = 
pi = a 1 =0. For this set of values of the parameters equation (|16[) admits the following explicit solutions 

y±(x) = l + C ± e ±2 ^ x . 

Let Y denote the set of all nonsymmetric solutions of p^|) such that the set of functions U introduced in Theorem lII.ll 
is linearly independent and S be the set of those elements in Y such that the potential is symmetric. The set S 
cannot be empty because if we take for instance k — 1/4, p = a = 0, pi = — a±, and y(x) = 1 + e x , then we obtain 
the following symmetric potential, namely 



(l + e*) 2 ' 

This shows that in general the potential Vj can be symmetric even though the coordinate transformation is not. 



III. NATANZON'S POTENTIAL 

In the same year as Iwata, Natanzon derived the most general potential such that the ODSE can be transformed 
into the hypergeometric equation. The starting point in his derivation is the following relation linking the potential 
V entering in the ODSE with the Bose invariant / 6j associated to the hypergeometric equation and the Schwarzian 
derivative {y,x} of the coordinate transformation [9J, namely 

(y) 2 I(y) + l{y,x} = k 2 -V(x), (43) 

with 

;|i , Jl-»i). 1 - a ) + M),^--l.,(l- 8 ) i mM*-*). (44, 

Here, the parameters /1, Ao and Ai are related to the parameters a, /?, and 7 of the hypergeometric equation as follows 

Ao = 7 — 1, Xi—a + f3~ 7, fi — fj — a. (45) 

Under the assumption that the coordinate transformation y does not depend on the energy of the particle and that 
the parameters /z, Ao and Ai are linear in fc 2 , i.e. 

1 - fj 2 = ak 2 - f, 1 - A| = c p k 2 -h p p = 0,l (46) 
Natanzon derived the following result 



9 



Theorem III.l The most general potential such that the ODSE can be transformed into the hypergeometric equation 
is given by 



V(x) = 



y 2 q-y) 2 

H 2 (y) 



o+(ci-oo)(2y-l) 5 A 



y(y- 1) 



4H(y) 



fy(y-l) + h {l-y) + h iy + l 

H(y) 



(47) 



with A = (a — Co — ci) 2 — 4cic and H(y) — ay 2 + (c! — c — a)j/ + co- Moreover, the coordinate transformation y = y(x) 
must satisfy the differential equation 



(y) 2 H(y) 
V(l -y) 2 



= 1. 



Proof. Let y = y{x) be the coordinate transformation leading from the hypergeometric equation to the ODSE. Such 
a transformation is determined by the condition (|4"5)l . By means of (l4"51) and (l4"4")l we obtain 



V(x) = k 



2 (y') 2 {k 2 [ay 2 + (c 1 -c -a)y + co]-[fy 2 + {h 1 -h - f)y + h }} 1 



-^{y,x}- 



(48) 



Assuming that (|4"rj]l holds and that y does not depend on the energy of the particle we end up with the following 
expression for the potential 

v ( x ) = k2 - ^ilyy ^ 2 W + ( c i - c ° - a )y + c °] - [fy 2 + fa - h o - f)y + h ]}- \{y, 

Let H(y) = ay 2 + (c\ — cq — a)y + cq. Then the above expression can be rewritten as 

The requirement that the potential does not depend on the energy of the particle implies that the coordinate trans- 
formation must satisfy the differential equation 



(y') 2 H(y) 
V(i - y) 2 

Hence, we obtain the following expression for the potential 



1. 



T// \ i; 1 , fy 2 + (fti - /to - f)y + feo 
y(x) = --{„, X } + _ . 



Using (|44p the above relation can be rewritten as 



4 \ y 



fy 2 + (hi - hp - f)y + h 
H{y) 



(49) 



Let us set g(y) — 4y 2 (l — y) 2 / H(y). Then, the ODE for y becomes (y ) 2 = g(y) and differentiating it we can rewrite 
the terms entering in the Schwarzian derivative as 

y" 



and thus 



y 2 dy 2 



4 U' 2y' Wg\dyJ 4dy 2 ' 



Hence, (PUR) can be rewritten as follows 



V(x) 3 ( dg 
16g \dy 



1 d 2 g fy 2 + (hi -h a - f)y + h 
Ady 2 H(y) 



(50) 
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Taking into account that 

dg_ = 8y(l - y)(l - 2y) Ay 2 (l-y) 2 
dy H(y) H 2 (y) 



G iy)i G(y) = 2ay + c% - cq - a, 



Jr/'V _ 64y 2 (l-y) 2 (l-2y) 2 64^(1 - y) 3 (l - 2y) | 16 ? /(1 - y) A ^ 



dy) H 2 (y) H 3 (y) yy ' H%y) 

:! (dg\ 2 _ 3(1 -2y) 2 3y(l - y)(l - 2y) c{y) | 3y 2 (l - yf 



8(6y 2 -6z/+l) 16y(l-y)(l-2y) 8a^ 2 (l - y) 2 8y 2 (l - j/) 2 2 

dy 2 H(y) m(y) U{V) H*{y) + H*{y) U [V) ' 



dgy ld 2 g _ 1 , y(l-y)(l-2y) c ^ | 2 a ?/ 2 (l-j / ) 2 5^ 2 (1 - Z/) 2 C 2 (y) 



we obtain 



16g\dyJ Adyi H(y) H 2 (y) m(y) AH 3 {y) 

Finally, substitution of the above expression into (f5T)]) gives 



V(x) 



H 2 (y) 



2a 



H(y) , (l-2y)G{y) 5G 2 (y) 



y 2 (l-y) 2 1/(1 -y) 4fr(y) 



/y 2 + (fci - fo - /)y + fop 

#(y) 



(51) 



By rewriting G 2 (y) in terms of A = (a — cp — ci) 2 — 4ciCp some trivial algebra leads to (|4"T)) . □ 



The next result shows the Natanzon class of potentials contains all classed of Iwata's potentials. This indicates that 
the method employed by Iwata is less general than that used by Natanzon. 

Theorem III. 2 Iwata's potentials \19\) , \32\) , and |ff7|) are a subclass of the Natanzon class U 



Proof. In Section [H] we derived formulae expressing the parameters a, b, and c in terms of the parameters entering 
in Iwata's potentials (TK))) . (13"2"1) . and (|3T|) . Substitution of these formulae into (|4l))) allows to express the parameters 
of Natanzon's potential in terms of those for the Iwata potentials. Hence, (|47p reproduces the potential Vi in (|T51) for 

Co = -, ho = -(a + cri), c x = 1 + P + ° , hi = -(1 + p + pi + a + 01), a = -, / = -1, (52) 



the potential V/j in (|32|) if 



co=~, ^o=0, ci = 1 + ^ + h 1 = -(p 2 + a 2 ), a= ^ f = ~P2, (53) 



and the potential Vjjj in (|37]) whenever 

Co = «o = -0-3, ci = - , /ii = -{p 3 + 03), a=-, f = -p 3 . (54) 

K K K 

At this point it is not difficult to construct a potential in the Natanzon class which does not belong to any Iwata class. 
For example, we can start by choosing choose / = 0, and ho = 1. This ensures that the potential cannot belong to 
the first or second Iwata class. Since we have hi = —pz — 03 = / + ho in ([54]) . we can choose h\ = 0, thus ensuring 
that this potential cannot be reproduced by any of the (flT)|) . (|3"2"|) . and (|3"T|) . Hence, Iwata's classes are contained in 
Natanzon's class. □ 



IV. GENERALIZED NATANZON POTENTIALS OF THE HEUN CLASS 

We shortly review some aspects of Milson's method [ll[ playing a fundamental role in the derivation of the most 
general potentials such that the ODSE can be reduced to a Heun equation or one of its confluent cases. Let us consider 
the second order linear equation 



. d 2 v , , . dv 
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By means of the Liouville transformation 

„(lA = pvn ( — / 

2a(t) 



we find that the standard form of (|55|) is 



d 2 / 



+ I(y)f(y) = 0, I(y) 



dy 2 ' " 4a 2 (y) 



4a(y)c(y) - 2a(y)£ + 2h(y)^ - 6 2 (y) 



(56) 



where J is the Bose invariant of ([55)1 (a, lilfl . If we make the coordinate transformation y = y(x) and set f(x) = f(y(x)), 
it can be seen that (|56|) becomes 

1 d2j y " dJ • mm = o, (57) 



(y') 2 dx 2 {y') 2 dx 

where the prime denotes differentiation with respect to x. The standard form of (|57|) can be attained by another 
Liouville transformation, namely 

J(x) = cxp ( - J ds \ u(x) = Vy 7 u{x) 



and we obtain 



dru i I 

— + J( x )u(x)=0, J(x) = (y) 2 I(y) + -{y,x}, (58) 



where the Schwarzian derivative is given by 

The solution of (|58l) will be expressed in terms of the solutions of (l55"j) by the following relation 

u(x) = (y )-V2 exp (^j ^ 

A comparison of the above expression with ((5|) show that the function s introduced by Iwata will be given by the 
following formula 

Equation (f58| will reduce to an ODSE whenever J(x) — k 2 — V(x). Hence, the potential V is entirely specified by 
the Bose invariant and the Schwarzian derivative of the coordinate transformation y. According to [ll| the potential 
V of the ODSE will not depend on k 2 if the Bose invariant admits a decomposition of the form I(y) — Ii(y)k 2 + Io(y) 
and the coordinate transformation y is a solution of the autonomous differential equation 

If this is the case, then the most general potential such that the ODSE can be transformed into (|55|) is 

v{x)= -m + — — ' (60) 

where a dot represents differentiation with respect to y. In what follows we extend the work of [ill ] by deriving 
expressions for the potential V in the case that (|55]l is a Heun equation, a confluent Heun equation, a double confluent 
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equation, a biconfluent Heun equation or a triconfluent Heun equation. We underline that Lemma 3 in [11| cannot 
be applied for the Heun equation and its confluent cases because it would require that in the Bose invariant 



T(y,k 2 ) 
A 2 (y,k 2 ) 



the polynomial A has has at most degree two, whereas the corresponding A for the Heun equation and its confluent 
cases has degree three. The next result shows that the most general potential associated to the Heun equation is 
controlled by ten parameters. 

Theorem IV. 1 Let /o, • • • , /t, go, ■ ■ ■ , <7 4 be real parameters. The most general potential such that the ODSE can be 
transformed into the Heun equation |3J) is given by 



V H (x) = 



y 2 (y~l) 2 (y-g) 2 

R 2 (y) 



R(y) 



G(y)R(y) - 2R(y)G(y) 

y(y - a) 



R(y)G 2 (v) 



5 (R(y)f 

y 2 (y-l) 2 (y-a) 2 4 R(y) 



S(y) 
R(y)' 



(61) 



where G(y) = 3y 2 — 2(a + l)y + a and 



R{y) = go{y - l) 2 (y - a) 2 + 9l y 2 (y - a) 2 + g 2 y 2 (y - I) 2 + g 3 y 2 {y - l)(y - a) + g 4 y{y - l)(y - a), 

S(y) = My - i) 2 (y - a) 2 + hy 2 (y - «) 2 + hy 2 {y - i) 2 + hy 2 (y - l)(y - a) + hv{y - l)(y - a). 

Moreover, the coordinate transformation y — y(x) satisfies the differential equation 

4y 2 (y-l) 2 (y-af 



(y 



R(y) 



Proof. By means of the relation e = a + f3 + l— 7 — 5 we can cast the Bose invariant of the Heun equation into the 
form 



i(y) = 



with 



My - l) 2 (y - a) 2 + W(y - a? + My 2 {y - l) 2 + hy 2 (y - i)(y - o) + Myjy - i)(y - o) 

% 2 (y- \) 2 {y- a) 2 
A = l-(l-7) 2 , A 1 = l-(1-J) 2 , A 2 = l-(a + /3- 7 -<5) 2 , 



(62) 



A 3 = 4a(3 - 2 7 c5 - 2(7 + 6) (a +0 + 1-^-5), A 4 = -Aq + 2aj5 + 2j(a + f3 + 1 - 7 - 8). 
Note that the Bose invariant is a rational function of the form 

77 \ = P ^ = c o + ci y + c ^y 2 + c zy 3 + C4 ^ 4 

{V) Q(y) Ay 2 (y-my-a) 2 

and it admits the decomposition I(y) = Ii(y)k 2 + Io(y) if and only if P(y) = R(y)k 2 + S(y) for some polynomials 
R(y) = a + a\y + a 2 y 2 + a 3 y 3 + a 4 y 4 and S(y) = b + b\y + b 2 y 2 + 63J/ 3 + b 4 y 4 . This will be the case if 



c = a k 2 + 6 = a 2 A , (63) 

ci = aik 2 + b x = -2a(a + 1)A + aA 4 , (64) 

c 2 = a 2 k 2 + b 2 = (a 2 + 4a + 1)A + a 2 Ai + A 2 + aA 3 - (a + 1)A 4 , (65) 

c 3 = a 3 k 2 + 63 = -2(a + 1)A - 2aAi - 2A 2 - (a + 1)A 3 + A 4 , (66) 

c 4 = a 4 fc 2 + 6 4 = A + A 4 + A 2 + A 3 + A 4 . (67) 



Note that (l63l implies that Ao is linear in k 2 and from (|64|) we conclude that A 4 is also linear in k 2 . We show now that 
Ai, A 2 and A3 are linear in k 2 . From ([65| -([67 | we see that a 2 Ai + A 2 + C1A3, — 2aAi — 2A 2 — (a+ 1)A3, and A 4 + A 2 + A3 
are also linear in k 2 . Let 

A1+A2+A3 = rik 2 + si, 
-2aAi - 2A 2 - (0+ 1)A 3 = r 2 k 2 + s 2 , 
a 2 X 1 +\ 2 + a\ 3 = r 3 k 2 + S3. 
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for some scalars n, • • ■ , ra, s±, • • • , S3. Let us rewrite the above non homogeneous system in matrix form, i.e. 

/11 1 \ f Xi \ f r A 

AA = k 2 T + B, A= \ -2a -1 -fa+1) , A = A 2 , T = r 2 , S = 







Since the parameter a in the Heun equation cannot assume the values and 1, we conclude that det(A) — (1 — a) 3 ^ 0. 
Hence, the matrix A is invertible and therefore Ai, A2, A3 are linear in k 2 . By setting A^ = gik 2 + fi with i = 0, 1, • • • ,4 
we can rewrite the Bose invariant (|62l) in the form I(y) = Ii(y)k 2 + Ia(y) with 



V(w-i) a (i/-o) 2 ' V(i/-i) 2 (i/-o) 3 ' 

and 

= 5o(y - i) 2 (y - a) 2 + <?iy 2 (y - a) 2 + 92y 2 (y - 1) 2 + g 3 y 2 (y - l)(y - a) + 54 y(y - i)(y - a), 

= fo(y - l) 2 (y - a) 2 + hy 2 (y - a) 2 + f 2 y 2 (y - I) 2 + hy 2 {y - \){y - a) + f 4 y(y - l)[y - a). 
Taking into account that 



we find that 



Moreover, 



(h) 2 . (R(y)) 2 R(y)G 2 (y) R{y)G{y) 



h 4y*(y-l) 2 (y-a) 2 R(y) y*(y - l) 4 (y - a) 4 y^y - l)*(y - a) 3 ' 



(68) 



f (,A = R{y) - R(y)G(y) + R(y)[3y-(a + l)} 3 R(y)G 2 (y) 

lW V(y_l)2( y _ a )2 y3 (y _ 1)3(y _ a) 3 + 2y\y~iy(y~ay [ ™> 

Finally, substituting (JBSJ) and (J55) into ([501) we obtain (JSTJ). □ 

In the next result we show that the Natanzon class is a special case of the potential fCTj) . 

Corollary IV. 2 Natanzon's class of solvable potentials |^7| ) is a subclass of the family of potentials h61\). 

Proof. The Bose invariant of the Heun equation will collapse into the Bose invariant of the hypergeometric equation 
whenever q = af3a and e = a + (3 + l — 7 — 6 = 0. In this case the parameters Ao, • • • , A4 in (1621) become 

A = l-(l- 7 ) 2 , Ai = l-(a + /3- 7 ) 2 , A 2 =0, A 3 = Aa(3 - 2 7 (a + (3 + 1 - 7), A 4 - -0A3. 

Since a denotes the position of a regular singular point in ([3]), we rename the parameters entering in Natanzon's 
potential according to the rule a — > a, Ci — > Ci, hi — > hi, and / — > /, where i = 0, 1. Then, in the notation of the 
previous theorem and making use of (|4"51) and P5|) we find that g^k 2 + f = A = Ao = c fc 2 — ho- Hence, it must be 
go = c and fo = — /lo- Moreover, <7ifc 2 + /1 = Ai = Ai = Cifc 2 — /ii and therefore we have 51 = c\ and /1 = — h\. 
Furthermore, A 2 = implies that gi = f% = 0. If we rewrite A3 as A3 = [1 — (a — /3) 2 ] + [1 — (a + (3— 7 ) 2 ] — [1 — (7— l) 2 ] 
and use (|4"5|) and (|4*51) we find that 

.93fc 2 + /3 = (a - co - cx)k 2 + h + hi - f. 

Hence, g% = a — oq — c\ and fs = ho + h\ — f. Finally, from the relation A4 = —0X3 we obtain g^ = —ag^ and 
fi = —a/3. This completes the proof. □ 

We consider now the case of the confluent Heun equation [l2[ 

d 2 v / 7 S \ dv Apay — a 



dy 2 \y y-l J dy y(y - 1) 
where 7, 5,p, a, a are arbitrary parameters. 



4 p t- + ( ^ <y) = °. ( 7 °) 
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Theorem IV. 3 Let o~q, ■ ■ ■ , 04, tq, ■ ■ ■ ,T4 be real parameters. The most general potential such that the ODSE can be 
transformed into the confluent Heun equation |7fl| ) is given by 



V C h{x) 



v 2 (y-i? 
R 2 (y) 



R(y) 



(2y-l)R(y) + iR(y) , (2y - l) 2 R{y) 5 (R(y)) 2 



2/(2/ — 1) 



y 2 (y~l) 2 4 J2(y) 



S(y) 
R(yY 



(71) 



where 



R(y) = <r Q y 2 + a 1 (y~l) 2 +a 2 y 2 (y-l) 2 + a 3 y 2 {y-l)+o- 4 y(y-l), 
S{y) = T y 2 + n{y-l) 2 + T 2 y 2 (y-1) 2 + T 3 y 2 {y-1)+T 4 y(y-1). 

Moreover, the coordinate transformation y — y(x) satisfies the differential equation 

■ 2 _Ay 2 (y-l) 2 



[y 



R{v) 



Proof. The Bose invariant of the confluent Heun equation can be written as 

A02, 2 + Xi(y - l) 2 + X 2 y 2 (y - 1) 2 + hy 2 (y - 1) + X 4 y( y - 1) + \r y( y - 1) 2 



4y 2 (y — l) 2 



(72) 



with 



A = 1 - (1 - 5) 2 , Ai = l-(l- 7 ) 2 , A 2 = -16p 2 ,A 3 = 8p(2a- 5), A 4 = -2(2<t + jS), X 5 = -8 P j. 
On the other side (j?2")) can also be expressed in terms of the Aj's as 



A 2 y 4 + (-2A 2 + A 3 + X 5 )y 3 + (A + Ai + A 2 - A 3 + A 4 - 2A 5 )y 2 + (-2Ai - A 4 + X 5 )y + X 1 



4y 2 {y-l) 2 

Moreover, I is a rational function of the form 

P(y) Co + + c 2 y 2 + c 3 y 3 + c 4 y 4 



Q(y) 



4y 2 (y - 1) 



and it admits the decomposition I(y) = Ii(y)k 2 + Io(y) if and only if P(y) = R(y)k 2 + S(y) for some polynomials 
R(y) = ao + aiy + a 2 y 2 + a 3 y 3 + a 4 y 4 and S(y) = b + b\y + b 2 y 2 + b 3 y 3 + b 4 y 4 . This will be the case if 

c = a k 2 + bo = Ai, 

ci = aifc 2 + 61 = — 2Ai — A 4 + A5, 

c 2 = a 2 k 2 + b 2 = Xo + Ai + A 2 — A 3 + X 4 — 2A 5 

C3 = a 3 k 2 + 63 = -2A 2 + A 3 + A 5 , 

c 4 = a 4 k 2 + 64 = A 2 . 

From the above set of equations it is straightforward to see that Ao, Ai, A 2 , A5 — A4, and A3 + A5 must be linear in 
k 2 . Let us rewrite the Bose invariant in terms of this quantities, that is 



i(y) 



X y 2 + X 1 (y- l) 2 + X 2 y 2 (y - l) 2 + (A 3 + X 5 )y 2 (y - 1) + (A 4 - X 5 )y(y - 1) 

4y 2 (2/-l) 2 



By setting A^ = cr^fc 2 + n with i = 0, 1,2 and A5 — A4 = a 3 k 2 + t 3i A3 + A5 = a 4 k 2 + t 4 we can rewrite (j?2")) in the 
form I(y) = I 4 (y)k 2 + Io(y) with 



h(y) 



R{y) 



■iy 2 (y- 1) 2 



. io(y) 



S(y) 



4yi(y-iy 



and 



R(y) = <r y 2 + <r 1 (y-l) 2 +<T 2 y 2 {y~l) 2 + <T 3 y 2 (y-l)+<T 4 y(y-l), 
S(y) = Toy 2 +n(y - l) 2 + T 2 y 2 (y - l) 2 + T 3 y 2 (y - 1) + T 4 y(y - 1). 
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Taking into account that 



we find that 



Moreover, 



h(y) 



R(y) (2y-l)R(y) 



iy 2 (y-l) 2 2y 3 (y-l) 3 ' 
(R(y)) 2 , (2y-l)*R(y) (2y - l)R{y) 



h 4y*(y - l) 2 R(y) y 4 (y - I) 4 y 3 (y - 1) 3 ' 
R(y) (2y-l)R(y) + R(y) , 3 (2y - l) 2 R(y) 



h[y) V(y-l) 2 y 3 (2/-l) 3 
Finally, substituting (J73J) and ([711) into ([101) we obtain ([7I]1. □ 
We consider now the case of the biconfluent Heun equation [l2| 



2 y 4 (y-i) 4 



(73) 



(74) 



d 2 i» / , N dt> 



( 7 -a-2)y- 



5 + (l + a)0 



v(y) = 0, 



(75) 



where a, /?, 7, 5 are arbitrary parameters. 



Theorem IV. 4 Let /3q, ■ ■ ■ , 0%, po, • • • , p% be real parameters. The most general potential such that the ODSE can be 
transformed into the biconfluent Heun equation |75|) is given by 



Vbch(x) 



y 



R 2 (y) 



R(y) + 



R(y) , R(y) 5 (R(y)) 2 



S(y) 
R(yY 



y y 2 4 R(y) 
where 

R{y) = fay 3 + foy 2 + Piy + Po, S{y) = -Ay 4 + p 3 y 3 + p 2 y 2 + fny + p . 
Moreover, the coordinate transformation y — y(x) satisfies the differential equation 

V 



(76) 



(y) 2 



R(vY 



Proof. The Bose invariant of the biconfluent Heun equation can be written as 

-Ay 4 - ABy 3 + (4 7 - B 2 )y 2 - 2Sy + 1 - a 2 



and it can be decomposed as I(y) = L 1 (y)k 2 + Io(y) if 



4y 2 



(77) 



B- -^k 2 - ^ 
P- 4 fc 4' 

It is straightforward to verify that 



A 1 -B 2 = (i 2 k 2 +p 2 , S = -^k 2 -^, I -a 2 = 6 k 2 +pa. 



R(y) S(y) 



4y 2 



with 



Considering that 



R{y) = fty 3 + foy 2 + Piy + B , S(y) = -Ay 4 + p 3 y 3 + p 2 y 2 + p x y + p a . 



ii(y) 



R(y) R(y) 

Ay 2 2y 3 '' 



we find that 



In addition, we have 



(A) 2 _ (R(y)) 2 | R(y) R(y) 

h iy 2 R(y) y 4 y 3 
R(y) R(y) , 3 R(y) 



4jr y 6 2 j/ 4 
Finally, substituting (JTHJ) and (HU) into flSO]) we obtain (JTHJ)- □ 
We consider now the case of the double confluent Heun equation [l2[ 



,d 2 v dv ( B_ 2 B-i 



a 2 '-y-J- + 1 + + B a + B iy + B 2 i/ v(y) = 0, 

«y dy \ y z y J 

where -B_2, ■ ■ • , Z?2 are arbitrary parameters. 

Theorem IV. 5 Let /iq, ■ ■ ■ , /X4, fo, • • • ,1/4 be real parameters. The most general potential such that the ODSE 
transformed into the double confluent Heun equation \8(fy is given by 



Vdch(x) 



y 



R 2 (y) 



R(y) 



2R(y) 5 (i?(y)) 5 



y 4 R(y) 



S(y) 
R(yY 



where 

R(y) = ^42/ 4 + M3y 3 + M22/ 2 + mJ + Mo, S(y) = v±y A + v 3 y 3 + v 2 y 2 + v x y + vq. 
Moreover, the coordinate transformation y — y[x) satisfies the differential equation 

(yf- V 



s(s) 



Proof. The Bose invariant of the double confluent Heun equation is 

4B 2 y 4 + 4B lV 3 + (AB + l)y 2 + 4B_ x y + 4B_ 2 



m = 

and it can be decomposed as I(y) = I\(y)k 2 + Io(y) if 



n M4 , 2 . v i D M3 , 2 . ^3 .d.i ; 2 , D MO , 2 , ^1 D MO , 2 . ^0 

B 2 = —k +—, B x = —k +— , 4B + l=/i 2 fc +i/ 2 , B_i = — fc + — , B_ 2 = — « 



It is straightforward to verify that 



with 



h(y) = ^, io(y)- S(y) 



-R(y) = M4y 4 + M3?/ 3 + M22/ 2 + fj-ty + Mo, ^(y) = ^ 4 y 4 + ^32/ 3 + ^y 2 + v x y + v . 

If we take into account that 

R{y) R(y) 



h{y) 



4y 4 y 



5 ' 



we arrive at 



We also have 



{h? (R(y)) 2 mv) g%) 

Zi 4y 4 i?(y) y6 y5 



h(y) 



R(y) 2R(y) | 5R(y) 

4y4 y5 y6 



Finally, substituting ([521) and into §U\) we obtain (|81[) . □ 
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We consider now the case of the triconfluent Heun equation [lj 



dy 2 

where Aq,Ai, A2 are arbitrary parameters. 



A + A lV + A 2 y 2 - ^y 4 ) ,<(,, \ = 0. 



(85) 



Theorem IV. 6 Let $i,uji with i = 0,1,2 be real parameters. The most general potential such that the ODSE can be 
transformed into the triconfluent Heun equation k85\) is given by 



Vtch(x) 



R 2 (y) 



4 R(y) 



S(y) 
R(yY 



(86) 



where 

R(y) = $2y 2 + #iv + i?o, S(y) = -9y 4 + uj 2 y 2 + my + w . 
Moreover, the coordinate transformation y — y(x) satisfies the differential equation 

(vf- 4 



R(y) 



Proof. The Bose invariant of the triconfluent Heun equation is trivially given by 

-9y 4 + AA 2 y 2 + AA iy + AA a 

i(y) = i ■ 

The decomposition I(y) = Ii(y)k 2 + Io(y) can be achieved if we choose 

A2 = tl k 2 + ^, Al = ^fe 2 + ^, A = %* + ( *>. 
4 4 ' 4 4 ' 4 4 



(87) 



Hence, we find that 



with 



Substituting 



into ([BO"]) we obtain 



h[y) = —r- , Io[y) = —r- 



R(y) = "92V 2 + $iy + #o, S(y) = -9y 4 + u 2 y 2 + my + wo- 



(A) 2 (R(y)) 2 T , , R(y) 

^T = 1^P h{y) = — 



□ 



We conclude this section by deriving the most general potential such that the ODSE can be transformed into the 
generalized Heun equation (|4]). Since this equation contains the Heun equation as a special case, the potential derived 
in the next result which depends upon 14 real parameters will also generalize the class of potentials obtained in 
Theorem UVU 

Theorem IV. 7 Let rji, ■ ■ ■ , tyr, £1, • • • ,£7 be real parameters. The most general potential such that the ODSE can be 
transformed into the generalized Heun equation is given by 



V G h(x) = 



y 2 (y- \) 2 {y-af 

R 2 (y) 



R(y) 



G(y)R{y) - 2R(y)G(y) 

y(y~ 



5 (R(y)) 2 

y 2 (y-l) 2 (y-a) 2 4 R{y) 



R(y)G 2 (y) 



S(y) 
R(yY 



where G(y) — 3y 2 — 2(a + l)y + a and 



R(y) = m(y - i) 2 (y - a) 2 + my 2 {y - a) 2 + my 2 {y - i) 2 + mv 2 (v - i) 2 (y - «) 2 
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+V5V(y - l) 2 (y - a) 2 + V6V 2 {y - 1)(V -a)+ r] 7 y(y - l)(y - a), 

S(y) = £1(2/ - i) 2 (y - «) 2 + 6y 2 (y - a) 2 + &y 2 (y - if + Uy 2 (y - i) 2 (y - «) 2 
- i) 2 (y - a) 2 + &y 2 (2/ - i)(y - a) + 62/(2/ - - a), 

Moreover, the coordinate transformation y = yix) satisfies the differential equation 

Kv ' R(y) 

Proof. Let us write the Bose invariant of the generalized Heun equation into the following form 
4y 2 (y-l) 2 (y-a) 2 I(y) = \ Q y(y~l)(y-a) + \ 1 y 2 (y-l)(y~a) + \ 2 y(y-l)(y~a) 2 + \ 3 y(y-l) 2 (y-a)+\ 4 y(y-l) 2 (y-a) 2 

+hy 2 (y - l)(y - «) 2 + W(y - l) 2 (y ~a) + \ 7 (y - l) 2 (y - a) 2 + X 8 y 2 (y - a) 2 + X g y 2 (y - I) 2 

X W y 2 (y - l) 2 (y - a) 2 + \ liy 3 (y - l)(y - a) (89) 

with 

\)=4/3 , Ai = 4/?i - 2ai« 2 , A 2 = -2a "i, A 3 = -2a a 2 , A 4 = -2a a 3 , A 5 = -2a 1 a 3 , 

\ e = —2a 2 a 3 , A7 = 2ao — a , A 8 = 2ai — a 2 , Ag = 2a 2 — a 2 , Ai = — a 2 , An = 4/? 2 . 
The Bose invariant is a rational function of the form 

P(y) c + c x y + c 2 y 2 + c 3 y 3 + c 4 y 4 + c 5 y 5 + c 6 y 6 



Q(y) Ay 2 {y-l) 2 {y-af 



and it admits the decomposition I(y) = Ii(y)k 2 + Io(y) if and only if P{y) = R{y)k 2 + S(y) for some polynomials 

R{y) = vo + my + my 2 + W + W + vsy 5 + my 6 and s{y) = £ + fry + i 2 y 2 + &y 3 + Uy 4 + &y 5 + &y 6 < This will 

be the case if 

Co = aok 2 + bo = a 2 \ 7 , c\ = a\k 2 + b\ = o\q — a 2 \ 2 — a\ 3 + a 2 A 4 — (2a + a 2 )A7, 
c 2 = a 2 k 2 +b 2 = -(l+a)\ Q +a\ 1 + (2a+a 2 )\ 2 + {l+2a)\ 3 -(2a+2a 2 )\ i -a 2 \ 5 -a\ 6 + (l+4a+a 2 )\ 7 +a 2 \$+\ 9 +a 2 \ w , 
c 3 = a 3 k 2 + 63 = A - (1 + a)Ai - (1 + 2a)A 2 - (2 + a)A 3 + (1 + 4a + a 2 )A 4 
+(2a + a 2 )A 5 + (1 + 2a)A 6 - (2 + 2a)A 7 - 2aA 8 - 2A 9 - (2a + 2a 2 )Ai + a\ u , 
c 4 = a 4 fc 2 + 6 4 = Ai + A 2 + A 3 - (2 + 2a)A 4 - (1 + 2a)A 5 - (2 + a)A 6 + A 7 + A 8 + A 9 + (1 + 4a + a 2 )Ai - (1 + a) An, 



c 5 = a 5 k 2 + 65 = A 4 + A 5 + A 6 — (2 + 2a)Ai + An, c 6 = a 6 k 2 + b 6 = A 
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It can be immediately seen that A7 and A10 are linear in k 2 . By applying the Gauss- Jordan elimination method to the 
remaining equations involving the Aj's we discover that Ag, Ag, A4 + A5 + A6 + A11, Ai + A 2 + A3 + A5 + aA6 + (1 + a)Au, 
and Ao — A3 — a(A 2 + A5 + Ag + An) must be also linear in k 2 . At this point the Bose invariant (I89[) can equivalently 
be written as 

4y 2 (y - l) 2 (y - a) 2 I(y) = X 7 (y ~ l) 2 (y - a) 2 + X 8 y 2 (y - a) 2 + X 9 y 2 (y - l) 2 + \ w y 2 {y - l) 2 (y - a) 2 



+ (A 4 + A 5 + A 6 + Xn)y(y - l) 2 (y - a) 2 + [Ai + A 2 + A 3 + A 5 + aA 6 + (1 + a)X u ]y 2 (y - l)(y - a) 
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[Ao - A 3 - a(A 2 + A 5 + A 6 + X n )]y(y - l)(y - a). 

Finally, if we set A 7 = r^k 2 + fi, A 8 = i] 2 k 2 + £2, Ag = r] 3 k 2 + f 3 , A10 = r\\k 2 + £4, A 4 + A 5 + A 6 + An = r] 5 k 2 + £5, 
Ai + A 2 + A 3 + A 5 + a\ e + (1 + a)An = i] 6 k 2 + £ 6 , and A - A3 - a(A 2 + A 5 + A 6 + An) = r^k 2 + £ 7 , where the % £,'s 
can be expressed in terms of the ci{, &,'s introduced before, we can rewrite the Bose invariant as I(y) — I\{y)k 2 + Io(y) 
where 

T I \ ?M T ( \ S ^ 



and 



Ay 2 {y-\) 2 {y-a) 2, Ki " 4y 2 (y - l) 2 {y - a) 2 ' 

R(y) = m(y - i) 2 (y - a) 2 + my 2 {y ~ a) 2 + my 2 (y - 1) 2 + mv 2 {y - i) 2 (y - «) 2 

+V5y(y - i) 2 (y - a) 2 + my 2 {y - i)(y - a) + ?77y(y - - a), 
s(y) = ti(v- !) 2 (y - «) 2 + &i/ a (i/ - «) 2 + &y 2 (y - !) 2 + ^y 2 (y - V 2 (y - «) 2 

H 5 y(y - i) 2 (y - a) 2 + &y 2 (y - i)(» - a) + bv(y - i)(» - a), 

At this point the derivation of the potential is the same as in Theorem lIV.ll □ 



V. CONCLUSIONS 



The Schrodinger equation is one of the key equations in Quantum Mechanics and plays a dominant role in many 
branches of science and technology. The history of analytical solvability of this equations has, of course, taken different 
routes. The concept of shape invariance is one example of a successful strategy [16]. Another line of attacking the 
problem is to transform a well known second order ordinary differential equation into the ODSE obtaining hereby a 
class of potentials. Such a procedure has been followed several times in the literature, mostly with the hypergeometric 
equation as a starting point. The Heun equation, its confluent forms and the generalized Heun equation offer the 
next possible generalization in this procedure. The class of Heun potentials allowing a transformation from the 
ODSE to one of the Heun equations has been derived in the paper. To appreciate how general these potentials 
are we mention that the potentials of Iwata class have five free parameters. The Natanzon class contains six free 
parameters. In comparison, the Heun class has as many as ten free parameters (with the exception of the biconflucnt 
Heun class which has eight parameters). The class of potentials of the generalized Heun equation has even fourteen 
free parameters. This fact will increase the number of possible applications. It is therefore worthwhile to extend the 
present study to consider the eigenvalue problem for bound states and other related problems. We intend to cover 
this and other topics in a forthcoming publication. 
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